The broken Alzheimer's disease genome.

The complex pathobiology of late-onset Alzheimer's disease (AD) poses significant challenges to therapeutic and preventative interventions. Despite these difficulties, genomics and related disciplines are allowing fundamental mechanistic insights to emerge with clarity, particularly with the introduction of high-resolution sequencing technologies. After all, the disrupted processes at the interface between DNA and gene expression, which we call the broken AD genome, offer detailed quantitative evidence unrestrained by preconceived notions about the disease. In addition to highlighting biological pathways beyond the classical pathology hallmarks, these advances have revitalized drug discovery efforts and are driving improvements in clinical tools. We review genetic, epigenomic, and gene expression findings related to AD pathogenesis and explore how their integration enables a better understanding of the multicellular imbalances contributing to this heterogeneous condition. The frontiers opening on the back of these research milestones promise a future of AD care that is both more personalized and predictive.

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants.

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Disregarding multimappers leads to biases in the functional assessment of NGS data.

BACKGROUND: Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. RESULTS: In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. CONCLUSION: Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Australian public perspectives on genomic newborn screening: which conditions should be included?

BACKGROUND: Implementing genomic sequencing into newborn screening programs allows for significant expansion in the number and scope of conditions detected. We sought to explore public preferences and perspectives on which conditions to include in genomic newborn screening (gNBS). METHODS: We recruited English-speaking members of the Australian public over 18 years of age, using social media, and invited them to participate in online focus groups. RESULTS: Seventy-five members of the public aged 23-72 participated in one of fifteen focus groups. Participants agreed that if prioritisation of conditions was necessary, childhood-onset conditions were more important to include than later-onset conditions. Despite the purpose of the focus groups being to elicit public preferences, participants wanted to defer to others, such as health professionals or those with a lived experience of each condition, to make decisions about which conditions to include. Many participants saw benefit in including conditions with no available treatment. Participants agreed that gNBS should be fully publicly funded. CONCLUSION: How many and which conditions are included in a gNBS program will be a complex decision requiring detailed assessment of benefits and costs alongside public and professional engagement. Our study provides support for implementing gNBS for treatable childhood-onset conditions.

Genomic analysis of a spontaneous unifoliate mutant reveals gene candidates associated with compound leaf development in Vigna unguiculata [L] Walp.

Molecular mechanisms which underpin compound leaf development in some legumes have been reported, but there is no previous study on the molecular genetic control of compound leaf formation in Vigna unguiculata (cowpea), an important dryland legume of African origin. In most studied species with compound leaves, class 1 KNOTTED-LIKE HOMEOBOX genes expressed in developing leaf primordia sustain morphogenetic activity, allowing leaf dissection and the development of leaflets. Other genes, such as, SINGLE LEAFLET1 in Medicago truncatula and Trifoliate in Solanum lycopersicum, are also implicated in regulating compound leaf patterning. To set the pace for an in-depth understanding of the genetics of compound leaf development in cowpea, we applied RNA-seq and whole genome shotgun sequence datasets of a spontaneous cowpea unifoliate mutant and its trifoliate wild-type cultivar to conduct comparative reference-based gene expression, de novo genome-wide isoform switch, and genome variant analyses between the two genotypes. Our results suggest that genomic variants upstream of LATE ELONGATED HYPOCOTYL and down-stream of REVEILLE4, BRASSINOSTERIOD INSENSITIVE1 and LATERAL ORGAN BOUNDARIES result in down-regulation of key components of cowpea circadian rhythm central oscillator and brassinosteroid signaling, resulting in unifoliate leaves and brassinosteroid-deficient-like phenotypes. We have stated hypotheses that will guide follow-up studies expected to provide more insights.

Genie: the first open-source ISO/IEC encoder for genomic data.

For the last two decades, the amount of genomic data produced by scientific and medical applications has been growing at a rapid pace. To enable software solutions that analyze, process, and transmit these data in an efficient and interoperable way, ISO and IEC released the first version of the compression standard MPEG-G in 2019. However, non-proprietary implementations of the standard are not openly available so far, limiting fair scientific assessment of the standard and, therefore, hindering its broad adoption. In this paper, we present Genie, to the best of our knowledge the first open-source encoder that compresses genomic data according to the MPEG-G standard. We demonstrate that Genie reaches state-of-the-art compression ratios while offering interoperability with any other standard-compliant decoder independent from its manufacturer. Finally, the ISO/IEC ecosystem ensures the long-term sustainability and decodability of the compressed data through the ISO/IEC-supported reference decoder.

OMD Curation Toolkit: a workflow for in-house curation of public omics datasets.

BACKGROUND: Major advances in sequencing technologies and the sharing of data and metadata in science have resulted in a wealth of publicly available datasets. However, working with and especially curating public omics datasets remains challenging despite these efforts. While a growing number of initiatives aim to re-use previous results, these present limitations that often lead to the need for further in-house curation and processing. RESULTS: Here, we present the Omics Dataset Curation Toolkit (OMD Curation Toolkit), a python3 package designed to accompany and guide the researcher during the curation process of metadata and fastq files of public omics datasets. This workflow provides a standardized framework with multiple capabilities (collection, control check, treatment and integration) to facilitate the arduous task of curating public sequencing data projects. While centered on the European Nucleotide Archive (ENA), the majority of the provided tools are generic and can be used to curate datasets from different sources. CONCLUSIONS: Thus, it offers valuable tools for the in-house curation previously needed to re-use public omics data. Due to its workflow structure and capabilities, it can be easily used and benefit investigators in developing novel omics meta-analyses based on sequencing data.

Pragmatic nationwide master observational trial based on genomic alterations in advanced solid tumors: KOrean Precision Medicine Networking Group Study of MOlecular profiling guided therapy based on genomic alterations in advanced Solid tumors (KOSMOS)-II study protocol KCSG AL-22-09.

BACKGROUND: Next-generation sequencing (NGS) has been introduced to many Korean institutions to support molecular diagnostics in cancer since 2017, when it became eligible for reimbursement by the National Health Insurance Service. However, the uptake of molecularly guided treatment (MGT) based on NGS results has been limited because of stringent regulations regarding prescriptions outside of approved indications, a lack of clinical trial opportunities, and limited access to molecular tumor boards (MTB) at most institutions. The KOSMOS-II study was designed to demonstrate the feasibility and effectiveness of MGT, informed by MTBs, using a nationwide precision medicine platform. METHODS: The KOSMOS-II trial is a large-scale nationwide master observational study. It involves a framework for screening patients with metastatic solid tumors for actionable genetic alterations based on local NGS testing. It recommends MGT through a remote and centralized MTB meeting held biweekly. MGT can include one of the following options: Tier 1, the therapeutic use of investigational drugs targeting genetic alterations such as ALK, EGFR, ERBB2, BRAF, FH, ROS1, and RET, or those with high tumor mutational burden; Tier 2, comprising drugs with approved indications or those permitted for treatment outside of the indications approved by the Health Insurance Review and Assessment Service of Korea; Tier 3, involving clinical trials matching the genetic alterations recommended by the MTB. Given the anticipated proportion of patients receiving MGT in the range of 50% ± 3.25%, this study aims to enroll 1,000 patients. Patients must have progressed to one or more lines of therapy and undergone NGS before enrollment. DISCUSSION: This pragmatic master protocol provides a mass-screening platform for rare genetic alterations and high-quality real-world data. Collateral clinical trials, translational studies, and clinico-genomic databases will contribute to generating evidence for drug repositioning and the development of new biomarkers. TRIAL REGISTRATION: NCT05525858.

Pancreatic ß-Cell Identity Change through the Lens of Single-Cell Omics Research.

The main hallmark in the development of both type 1 and type 2 diabetes is a decline in functional ß-cell mass. This decline is predominantly attributed to ß-cell death, although recent findings suggest that the loss of ß-cell identity may also contribute to ß-cell dysfunction. This phenomenon is characterized by a reduced expression of key markers associated with ß-cell identity. This review delves into the insights gained from single-cell omics research specifically focused on ß-cell identity. It highlights how single-cell omics based studies have uncovered an unexpected level of heterogeneity among ß-cells and have facilitated the identification of distinct ß-cell subpopulations through the discovery of cell surface markers, transcriptional regulators, the upregulation of stress-related genes, and alterations in chromatin activity. Furthermore, specific subsets of ß-cells have been identified in diabetes, such as displaying an immature, dedifferentiated gene signature, expressing significantly lower insulin mRNA levels, and expressing increased ß-cell precursor markers. Additionally, single-cell omics has increased insight into the detrimental effects of diabetes-associated conditions, including endoplasmic reticulum stress, oxidative stress, and inflammation, on ß-cell identity. Lastly, this review outlines the factors that may influence the identification of ß-cell subpopulations when designing and performing a single-cell omics experiment.

Genomic insights into the diversity, virulence, and antimicrobial resistance of group B Streptococcus clinical isolates from Saudi Arabia.

Introduction: Detailed assessment of the population structure of group B Streptococcus (GBS) among adults is still lacking in Saudi Arabia. Here we characterized a representative collection of isolates from colonized and infected adults. Methods: GBS isolates (n=89) were sequenced by Illumina and screened for virulence and antimicrobial resistance determinants. Genetic diversity was assessed by single nucleotide polymorphisms and core-genome MLST analyses. Results: Genome sequences revealed 28 sequence types (STs) and nine distinct serotypes, including uncommon serotypes VII and VIII. Majority of these STs (n=76) belonged to the human-associated clonal complexes (CCs) CC1 (33.71%), CC19 (25.84%), CC17 (11.24%), CC10/CC12 (7.87%), and CC452 (6.74%). Major CCs exhibited intra-lineage serotype diversity, except for the hypervirulent CC17, which exclusively expressed serotype III. Virulence profiling revealed that nearly all isolates (94.38%) carried at least one of the four alpha family protein genes (i.e., alphaC, alp1, alp2/3, and rib), and 92.13% expressed one of the two serine-rich repeat surface proteins Srr1 or Srr2. In addition, most isolates harbored the pilus island (PI)-2a alone (15.73%) or in combination with PI-1 (62.92%), and those carrying PI-2b alone (10.11%) belonged to CC17. Phylogenetic analysis grouped the sequenced isolates according to CCs and further subdivided them along with their serotypes. Overall, isolates across all CC1 phylogenetic clusters expressed Srr1 and carried the PI-1 and PI-2a loci, but differed in genes encoding the alpha-like proteins. CC19 clusters were dominated by the III/rib/srr1/PI-1+PI-2a (43.48%, 10/23) and V/alp1/srr1/PI-1+PI-2a (34.78%, 8/23) lineages, whereas most CC17 isolates (90%, 9/10) had the same III/rib/srr2/P1-2b genetic background. Interestingly, genes encoding the CC17-specific adhesins HvgA and Srr2 were detected in phylogenetically distant isolates belonging to ST1212, suggesting that other highly virulent strains might be circulating within the species. Resistance to macrolides and/or lincosamides across all major CCs (n=48) was associated with the acquisition of erm(B) (62.5%, 30/48), erm(A) (27.1%, 13/48), lsa(C) (8.3%, 4/48), and mef(A) (2.1%, 1/48) genes, whereas resistance to tetracycline was mainly mediated by presence of tet(M) (64.18%, 43/67) and tet(O) (20.9%, 14/67) alone or in combination (13.43%, 9/67). Discussion: These findings underscore the necessity for more rigorous characterization of GBS isolates causing infections.

Estimating genomic relationships of metafounders across and within breeds using maximum likelihood, pseudo-expectation-maximization maximum likelihood and increase of relationships.

BACKGROUND: The theory of "metafounders" proposes a unified framework for relationships across base populations within breeds (e.g. unknown parent groups), and base populations across breeds (crosses) together with a sensible compatibility with genomic relationships. Considering metafounders might be advantageous in pedigree best linear unbiased prediction (BLUP) or single-step genomic BLUP. Existing methods to estimate relationships across metafounders Γ are not well adapted to highly unbalanced data, genotyped individuals far from base populations, or many unknown parent groups (within breed per year of birth). METHODS: We derive likelihood methods to estimate Γ . For a single metafounder, summary statistics of pedigree and genomic relationships allow deriving a cubic equation with the real root being the maximum likelihood (ML) estimate of Γ . This equation is tested with Lacaune sheep data. For several metafounders, we split the first derivative of the complete likelihood in a term related to Γ , and a second term related to Mendelian sampling variances. Approximating the first derivative by its first term results in a pseudo-EM algorithm that iteratively updates the estimate of Γ by the corresponding block of the H-matrix. The method extends to complex situations with groups defined by year of birth, modelling the increase of Γ using estimates of the rate of increase of inbreeding ( Δ F ), resulting in an expanded Γ and in a pseudo-EM+ Δ F algorithm. We compare these methods with the generalized least squares (GLS) method using simulated data: complex crosses of two breeds in equal or unsymmetrical proportions; and in two breeds, with 10 groups per year of birth within breed. We simulate genotyping in all generations or in the last ones. RESULTS: For a single metafounder, the ML estimates of the Lacaune data corresponded to the maximum. For simulated data, when genotypes were spread across all generations, both GLS and pseudo-EM(+ Δ F ) methods were accurate. With genotypes only available in the most recent generations, the GLS method was biased, whereas the pseudo-EM(+ Δ F ) approach yielded more accurate and unbiased estimates. CONCLUSIONS: We derived ML, pseudo-EM and pseudo-EM+ Δ F methods to estimate Γ in many realistic settings. Estimates are accurate in real and simulated data and have a low computational cost.

Redefining and interpreting genomic relationships of metafounders.

Metafounders are a useful concept to characterize relationships within and across populations, and to help genetic evaluations because they help modelling the means and variances of unknown base population animals. Current definitions of metafounder relationships are sensitive to the choice of reference alleles and have not been compared to their counterparts in population genetics-namely, heterozygosities, FST coefficients, and genetic distances. We redefine the relationships across populations with an arbitrary base of a maximum heterozygosity population in Hardy-Weinberg equilibrium. Then, the relationship between or within populations is a cross-product of the form Γ b , b ' = 2 n 2 p b - 1 2 p b ' - 1 ' with p being vectors of allele frequencies at n markers in populations b and b ' . This is simply the genomic relationship of two pseudo-individuals whose genotypes are equal to twice the allele frequencies. We also show that this coding is invariant to the choice of reference alleles. In addition, standard population genetics metrics (inbreeding coefficients of various forms; FST differentiation coefficients; segregation variance; and Nei's genetic distance) can be obtained from elements of matrix Γ .

Christensenella strain resources, genomic/metabolomic profiling, and association with host at species level.

The gut commensal bacteria Christensenellaceae species are negatively associated with many metabolic diseases, and have been seen as promising next-generation probiotics. However, the cultured Christensenellaceae strain resources were limited, and their beneficial mechanisms for improving metabolic diseases have yet to be explored. In this study, we developed a method that enabled the enrichment and cultivation of Christensenellaceae strains from fecal samples. Using this method, a collection of Christensenellaceae Gut Microbial Biobank (ChrisGMB) was established, composed of 87 strains and genomes that represent 14 species of 8 genera. Seven species were first described and the cultured Christensenellaceae resources have been significantly expanded at species and strain levels. Christensenella strains exerted different abilities in utilization of various complex polysaccharides and other carbon sources, exhibited host-adaptation capabilities such as acid tolerance and bile tolerance, produced a wide range of volatile probiotic metabolites and secondary bile acids. Cohort analyses demonstrated that Christensenellaceae and Christensenella were prevalent in various cohorts and the abundances were significantly reduced in T2D and OB cohorts. At species level, Christensenellaceae showed different changes among healthy and disease cohorts. C. faecalis, F. tenuis, L. tenuis, and Guo. tenuis significantly reduced in all the metabolic disease cohorts. The relative abundances of C. minuta, C. hongkongensis and C. massiliensis showed no significant change in NAFLD and ACVD. and C. tenuis and C. acetigenes showed no significant change in ACVD, and Q. tenuis and Geh. tenuis showed no significant change in NAFLD, when compared with the HC cohort. So far as we know, this is the largest collection of cultured resource and first exploration of Christensenellaceae prevalences and abundances at species level.

The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans.

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

Genomic landscape and distinct molecular subtypes of primary testicular lymphoma.

Primary testicular lymphoma (PTL) is a rare lymphoma predominantly occurring in the elderly male population. It is characterized by a limited response to treatment and a heightened tendency towards relapse. Histologically, approximately 90% of PTL cases are classified as diffuse large B-cell lymphomas (DLBCL). Genetic features of PTL were delineated in a limited scope within several independent studies. Some of the articles which analyzed the genetic characterization of DLBCL have incorporated PTL samples, but these have been constrained by small sample sizes. In addition, there have been an absence of independent molecular typing studies of PTL. This report summarizes the common mutational features, copy number variations (CNVs) and molecular typing of PTL patients, based on whole-exome sequencing (WES) conducted on a cohort of 25 PTL patients. Among them, HLA, CDKN2A and MYD88 had a high mutation frequency. In addition, we found two core mutational characteristics in PTL including mutation in genes linked to genomic instability (TP53 and CDKN2A) and mutation in immune-related genes (HLA, MYD88, CD79B). We performed molecular typing of 25 PTL patients into C1 subtype with predominantly TP53 mutations and C2 subtype with predominantly HLA mutations. Notably, mutations in the TP53 gene predicted a poor outcome in most types of lymphomas. However, the C1 subtype, dominated by TP53 mutations, had a better prognosis compared to the C2 subtype in PTL. C2 subtype exhibited a worse prognosis, aligning with our finding that the mechanism of immune escape in PTL was primarily the deletions of HLA rather than PD-L1/PD-L2 alterations, a contrast to other DLBCLs. Moreover, we calculated the tumor mutation burden (TMB) and identified that TMB can predict prognosis and recurrence rate in PTL. Our study underscores the significance of molecular typing in PTL based on mutational characteristics, which plays a crucial role in prognostication and guiding therapeutic strategies for patients.

Comparative genomic analyses provide insight into the pathogenicity of three Pseudomonas syringae pv. actinidiae strains from Anhui Province, China.

BACKGROUND: Pseudomonas syringae pv. actinidiae (Psa) is an important bacterial plant pathogen that causes severe damage to the kiwifruit industry worldwide. Three Psa strains were recently obtained from different kiwifruit orchards in Anhui Province, China. The present study mainly focused on the variations in virulence and genome characteristics of these strains based on the pathogenicity assays and comparative genomic analyses. RESULTS: Three strains were identified as biovar 3 (Psa3), along with strain QSY6 showing higher virulence than JZY2 and YXH1 in pathogenicity assays. The whole genome assembly revealed that each of the three strains had a circular chromosome and a complete plasmid. The chromosome sizes ranged from 6.5 to 6.6 Mb with a GC content of approximately 58.39 to 58.46%, and a predicted number of protein-coding sequences ranging from 5,884 to 6,019. The three strains clustered tightly with 8 Psa3 reference strains in terms of average nucleotide identity (ANI), whole-genome-based phylogenetic analysis, and pangenome analysis, while they were evolutionarily distinct from other biovars (Psa1 and Psa5). Variations were observed in the repertoire of effectors of the type III secretion system among all 15 strains. Moreover, synteny analysis of the three sequenced strains revealed eight genomic regions containing 308 genes exclusively present in the highly virulent strain QSY6. Further investigation of these genes showed that 16 virulence-related genes highlight several key factors, such as effector delivery systems (type III secretion systems) and adherence (type IV pilus), which might be crucial for the virulence of QSY6. CONCLUSION: Three Psa strains were identified and showed variant virulence in kiwifruit plant. Complete genome sequences and comparative genomic analyses further provided a theoretical basis for the potential pathogenic factors responsible for kiwifruit bacterial canker.

AnnoView enables large-scale analysis, comparison, and visualization of microbial gene neighborhoods.

The analysis and comparison of gene neighborhoods is a powerful approach for exploring microbial genome structure, function, and evolution. Although numerous tools exist for genome visualization and comparison, genome exploration across large genomic databases or user-generated datasets remains a challenge. Here, we introduce AnnoView, a web server designed for interactive exploration of gene neighborhoods across the bacterial and archaeal tree of life. Our server offers users the ability to identify, compare, and visualize gene neighborhoods of interest from 30 238 bacterial genomes and 1672 archaeal genomes, through integration with the comprehensive Genome Taxonomy Database and AnnoTree databases. Identified gene neighborhoods can be visualized using pre-computed functional annotations from different sources such as KEGG, Pfam and TIGRFAM, or clustered based on similarity. Alternatively, users can upload and explore their own custom genomic datasets in GBK, GFF or CSV format, or use AnnoView as a genome browser for relatively small genomes (e.g. viruses and plasmids). Ultimately, we anticipate that AnnoView will catalyze biological discovery by enabling user-friendly search, comparison, and visualization of genomic data. AnnoView is available at http://annoview.uwaterloo.ca.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.